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ABSTRACT 

The inner degrees of the Galactic center contain a large population of filamentary structures observed 
at radio frequencies. These so-called non-thermal radio filaments (NRFs) trace magnetic field lines 
and have attracted significant interest due to their hard (S„ oc v ~ 01 ± a4 ) synchrotron emission 
spectra. The origin of these filaments remains poorly understood. We show that the electrons and 
positrons created through the annihilations of a relatively light (~5-10 GeV) dark matter particle 
with the cross section predicted for a simple thermal relic can provide a compelling match to the 
intensity, spectral shape, and flux variation of the NRFs. Furthermore, the characteristics of the dark 
matter particle necessary to explain the synchrotron emission from the NRFs is consistent with those 
required to explain the excess 7-ray emission observed from the Galactic center by the Fermi-LAT, as 
well as the direct detection signals observed by CoGeNT and DAMA/LIBRA. 



1. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) pro- 
vide an attractive class of candidates for the dark matter 



of our universe (Bertone et al. 2005 Bergstrom 2009 1 . 
The WIMP paradigm is motivated in part by the re- 
alization that particles with weak-scale interactions and 
masses will naturally freeze-out in the early universe with 
a relic a bundance similar to the observed density of dark 

an observation referred to as 



matter (Steigman 
the "WIMP Mirac 



1979], 

Barring any complicating factors 



such as conannihilations, resonances, or S'-wave suppres 
sion, dark matter candidates motivated by the WIMP 
miracle annihilate with a cross section of approximately 
(TV ~ 3 x 10~ 26 cm 3 s _1 (where v is the relative velocity of 
the annihilating WIMPs). These annihilations produce 
stable particles, including photons at 7-ray energies, as 
well as protons, electrons, neutrinos, and their antipar- 
ticles. Because the dark matter annihilation rate scales 
with the square of the dark matter density, regions such 
as the Galactic center, dwarf spheroidal galaxies, and 
galaxy clusters represent promising locations for searches 
for indirect signatures of dark matter. 

The launch of the Fermi Large Area Telescope (Fermi- 
LAT) in 2008 has greatly ex panded our view of the 7- 
ray sky (Atwood et al. 20091. In addition to its signif- 
icantly enhanced effective area, the unparalleled spatial 
and energy resolution of the Fermi-LAT has allowed for 
the separation of point sources in the Galactic center, 
revealing much more detailed informa tion about the dif- 
fuse 7-ray em ission from this r egion ([A hdoet^aL 20091 
V it ale et al. [2009 ) . Recently, |Hooper fc Uoodenough 
( |201ip identified an excess of 7-rays within approxi- 
mately 175 pc of the Galactic center in the energy range 
of ^500 MeV to 7 GeV, and showed that this could be 
explained by 7-10 GeV dark matter particles annihilat- 
ing into t + t~ pairs, possibly among other leptonic final 
states. This range of dark matter masses also provides a 
suitable match to the direct detection signals observed by 
the DAMA/LIBRA ( |Bernabei et al.||2010| |Hooper et ah 



2010 [) and CoGeNT collaborations flAalseth et al.J 2011 



Aalscth et al. 2010 



Hooper fc Kel1o1[2uTT1). Whi e up 



coming limits from CMB constraints (platyer e t al. 12009 



Hutsi et al.||2011 ) as well as LEP constraints ( jl'ox et a 
201 1| ), error bars in these measurements, along with un- 



certainties in the dark matter density both locally and 
near the galactic center, are more than sufficient to rem- 
edy any tension between these models. 

If annihilating dark matter particles are in fact respon- 
sible for the flux of 7-rays observed from the Galactic 
center, this signal is guaranteed to be accompanied by 
the production of a hard population of electrons and 
positrons carrying a significant percentage of the total 
annihilation energy. In the case of democratic annihi- 
lation (equal number of annihilations to each family of 
charged leptons: e + e~, and t + t~), the combined 

electrons and positrons carry away nearly an order of 
magnitude more energy than the 7-rays observable by 
the Fermi-LAT (8.4 GeV per annihilation, as compared 
to 1.1 GeV into 7-rays for a 8 GeV WIMP). For this 
reason, the synchrotron radiation from the electrons and 
positrons produced in dark matter annihilations provides 
a particularly promising test of the dark matter interpre- 
tation of the 7-ray flux from the Galactic center. Along 
these lines, it has been shown that the characteristics 
of the galactic synchrotron excess known as the WMAP 
Haze are consistent with scenarios in which annihilating 
dark matt er is the source of the ob served Galactic cen- 
ter 7-rays ( Hooper fc Lindeii]|2011 1. 



In addition to the hig h dark matter density expected 
in the Galactic center (Navarro et al. 19961) , galactic 



studies of dark matter synchrotron are enticing due to 
the many survey observations und ertaken over severa l 
decades. Observations at 74 MHz ( Brogan et al.| 2003|), 
330 MHz (fPedlar et al.|1989[|Ananthar arnaiah et a lTT9"9" l 



|LaRosa et al.||2000| |Nord et al.|[2004|, 1.4 G Hz ( |Yusef 



|/adeh et al. 20 0l|~5 GHz fEaRosa et al. 2004) and higher 
frequencies ( jSofue et al.|[l9 86; Reich et al. 2000) allow 
for the additional modeling of spectral features, in ad- 



2 



dition to the spatial characteristics, of radio sources. A 
particularly interesting probe for synchrotron signals of 
dark matter annihilation concerns the population of long 
(—40 pc) and thin (~1 pc) filamentary structures located 
between 10 and 200 pc from the Galactic center, which 



et al. 
radio 



have been i dentified at radio wavelengths (|Yusef-Zadeh 

T 



1984). These structures, known as non-thermal 
laments (NRFs), are characterized primarily by 
their extremely hard radio spectra, with a spectral index 
a between -0.5 and +0.3 (where S„ ~ ^ Q ), as well as 
their preferential alignment perpendicular to the Galac- 
tic Plane. In several notable cases such as the Radio Arc, 
high resolution surveys have found a tangled network of 
separate NRFs w hich contribute to the ove rall emission 
in the structure (Yusef-Zadeh et al.l 1984). The radio 



emission from these filaments is highly polarized, imply- 
ing that the radio sources are powered by synchrotron 
emission in a highl y ordered poloidal magnetic field of 
strength -100 fiG (Tsuboi et aL]p86 |. 



While several dynamical mechanisms have been pro- 
posed to explain the comparatively strong magnetic fields 



netic field 


Rosner & Bodo 


L996 


Boldyrev & Yusef- 


Zadeh 


200G 


, the hard synchrotron spectrum observed 



served synchrotron spectral index, a, stemming from a 
power-law electron injection spectrum with index p can 
be modeled as p = 2a — 1, where p is the power law in- 
dex of the electron injection spectrum and a is the power 
law index of the output synchrotron spectrum. Thus, the 



shock acceleration (Blandford & Eichler||1987 


Jones & 


Ellison 


1991 


Malkov fe O'C Drury||2001|) would be ex- 



than a ~ -0.7, which is significantly softer than obser- 
vations of NR Fs such as GO. 2-0.0 , commonly known as 
the Radio Arc. Lesch et al. ( 1988 ) modeled the observed 
synchrotron spectrum of the Radio Arc, finding the syn- 
chrotron spectrum to be best fit by an essentially mo- 
noenergetic electron spectrum at approximately 7 GeV. 
Further observations have found a ~10 GHz turnover 
in the synchrotron spectrum of many NRFs, implying 
a strongly peaked electron energy spectrum at several 
GeV pr opagating in a magnetic field on the order of 
100 juG ( IBoldyrev fe Yusef- Zadeh|[2006l ). 

Several mechanisms have been proposed to explain 
the origin of such a hi ghly peaked electron s pectrum in 
NRFs. Mo st not ably, |Lesch fe Reich| ( [T992| , and later 
advocated a pile-up scenario in which 



Lieb et al. 



( |2004D , 



magnetic reconnection zones formed in collisions between 
NRFs and molecular clouds create electrical potentials 
which can accelerate electrons to a single energy speci- 
fied by the relative strength of the electromagnetic po- 
tential and the synchrotron energy loss rate. This ex- 
planation, however, is problematic for three reasons: (1) 
not all bri ght NRFs have observed molecular cloud as - 
sociations ( Lang et al.|1999a Law fe Yusef-Zadeh|2 004), 
(2) recent simulations of coilisional reconnection regimes 
imply a maximum electron energy of less than 10 McV, 
several orders of magnitude below th at needed to ex- 



plain the observed synchr otron signal (Lyubarsky 2005 
Zanotti & Dumbser 2011), although collisionless recon- 
nection, already proposed in the case of 7-ray bursts, 
may be significantly more effective in accelerating leptons 



to hig h energies ( McKinney fe Uzdensky|[2010| Lazarian 
et al. 2011), and (3) the maximum energy 01 particles 
moving through magnetic reconnection is expected to 
depend sensitively on the ge ometry of the reconnection 
region (Lazarian et al. 2011), which would be expected 
to produce very different synchrotron emission spectra 
in different NRFs. An alternate scenario involves elec- 



tron acceleration in star formation regions (Yusef-Zadeh 



20031, although the necessary mechanism to create a 
peaked electron spectrum in star formation regions is not 
known. 

In this paper, we propose that dark matter annihi- 
lation may provide a physical basis for the nearly mo- 
noenergetic electron spectrum necessary to explain NRF 
emission. Specifically, we find that the highly peaked 
electron injection spectrum naturally produced by an- 
nihilations of light dark matter particles correctly pro- 
duces the hard and bright synchrotron emission spec- 
trum observed in multiple NRFs. In § [2j we review the 
astrophysics of Alimentary arcs and demonstrate that the 
ordered magnetic field structure within these objects re- 
quires that their synchrotron emission results from elec- 
trons and positrons injected from within the filaments 
themselves, rather than from external sources. In § [3j 
we show that relatively light dark matter parti cles an- 
nihilating to leptons, such as those proposed by |Hooper] 
& Goodenough (2011) to explain the Galactic center 7- 
ray excess, are also predicted to inject electrons with a 
spectrum and intensity naturally capable of explaining 
the synchrotron intensity and spectrum observed from 
NRFs. In § |4j we show that synchrotron emission pro- 
duced as a product of dark matter annihilation can ex- 
plain the characteristics of the NRF population, includ- 
ing their spectral conformity and the spherical symmetry 
of their intensity with respect to the Galactic center. Fi- 
nally, in § [5j we present several testable predictions for 
this dark matter annihilation scenario with the goal of 
discriminating this possibility from astrophysical source 
mechanisms. Further study of the Milky Way's radio 
filaments may play a critical role in untangling various 
interpretations of the 7-ray excess observed in the Galac- 
tic center. 

2. THE ASTROPHYSICS OF THE NON-THERMAL RADIO 
FILAMENTS 

The strength of the magnetic fields in NRFs have 
been estimated through a variety of means. Early es- 
timates centered on the brightest NRF, GO. 2-0.0 (the 
Radio Arc), and were based on a comparison between 
the magnetic field pressure and the estimated ram pres- 
sure from nearby molecular cloud interactions, w hich in- 
dicated magnetic field strengths as high as 1 mG ( Yusef- 



Zadeh fe Morris|198 7 
recent observations 



Morris fe Yusef-Zadeh 1989[ ) I More 
lowever, have pointed to somewhat 



weaker magnetic fields (~ 100 /j,G) among the popula- 
tion of NRFs for three reasons: (1) the observations of 
kinks in several NRFs such as G359. 1-0.2 (the Snake) 
and G358.85+0.47 (the Pelican) imply that the mag- 
netic field pres sure may be significantly smaller than the 
ram pressure (Gray et al. 1995[ Lang et al. 1999a), (2) 
the compact radial extent of the filaments is difficult 
to explain if a 1 mG magnetic field is surrounded by 
a reg ion with a significa ntly weaker galactic magnetic 
field ( |LaRosa et al.|2004| , and (3) synchrotron models of 



3 



the radio spectrum imply equipartition magnetic fields 



between approx imately 50 and 200 /iG ( Anantharama- 
iah et al. 1991). Furthermore, in the overall population 



of JNKFs, a turnover of the hard synchrotron spectrum at 
~10 GHz is observed, i mplying; a magnetic field strength 
on the order of 100 /iG ( |Boldyrev fc Yusef-Zadehff 2006 ) . 

Although it has been suggested that these relatively 
strong magnetic fields may simply trace an extremely 
strong poloidal magnetic field formed in the Galactic 
center during the proto-halo phase of galaxy formation 
flSofue fc Fujimoto||1987| |Chandran et al.||2000| [Crocker 



et al.|2UlTi ), more recent studies have instead interpreted 



JNKFs to contain local enha ncements of the relatively 



weak diffuse magne tic field ( Boldyrev fc Yusef-Zadeh 
20061 iFerrierel 120091 and references therein.). ILaRosa 



et al. Q2005 ) studied the synchrotron spectral index and 
flux densities at 74 and 330 MHz, and calculated the 
strength of the large-scale diffuse magnetic field to be 
approximately ~10 /iG in the inner Milky Way. Further- 
more, a ubiquitous poloidal field is unable to explain re- 
cent observations of several NRFs w hich are not aligned 



perpendicular to the Galactic Disk (Lang et al. 1999a 



Nord et al.|2004 ). Finally, depolarization due to Faraday 



rotation also implies the strength of the diffu se magnetic 
field along the line-of-sight to be ~7 /iG ( |Gray et al 



19951. 

The strong magnetic fields in NRFs relative to those in 
the surrounding galactic medium can greatly affect the 
propagation of cosmic rays into and out of the filaments. 
The high degree of polarization of the synchrotron emis- 
sion from NRF s implies that their magnet ic fields are 

Assuming 



highly ordered (Yusef-Zadeh & Morris 



1987) 



that depolarization is dominated by turbulence in the 
magnetic field, the observed fractional polarization can 
be approximated as: 



P ~Po(l) 



Bl 



pod) 



4 

Brp 



(1) 



where po(7) is the intrinsic polarization (determined by 
the index 7 of the electron injection spectrum), Bq is the 
intensity of the ordered field, B r is the intensity of a ran- 

and B T = y/~B$ 



Bl 



is 



dom field with variance |S ; 

the total field. This solution stands as an equality in the 
case of a synchrotron spect rum a = -1 , and varies only 
slightly in all relevant cases ( Burn|1966 1 . Observations of 
~60% polarization in multiple JNKFs ([Yusef-Zadeh et al7 
|2004[ ) thus indicate magnetic fields with ordered tractions 
(Bq/Bj,) greater than ^80%. In reality, depolarizing ef- 
fects such as Fa raday rotation may dominate synchrotron 
depolarization ( Lang et al.||1999a ), implying a magnetic 
field which is significantly more ordered than required by 
this lower limit. 

In the presence of highly ordered magnetic fields, the 
gyroradius of electrons impinging perpendicular to the 
field lines of the NRFs is extremely small, leading the 
filaments to act as a magnetic mirror which quickly re- 
jects incident electrons. In Fig. [TJ we show the aver- 
age time that an electron impacting a cylindrical mag- 
netic field with total strength 100 /iG remains within a 
filament before being evicted, as a function of the or- 
dered fraction of the magnetic field. We employ a fila- 
ment radius of 0.412 pc (10), and imagine the filament 
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Fig. 1. — The average time that 10 GeV electrons spend within a 
filament of total magnetic field of 100 fiG after impacting its side, 
as a function of the ordered fraction of the magnetic field energy 
density within the filament. This time is insignificant compared 
to the ~10 12 second energy loss time for these particles, indicat- 
ing that externally produced electrons contribute only a negligible 
fraction of the synchrotron radiation from NRFs. 

as an infinitely long tube in order to safely ignore any 
edge effects. Using Monte Carlo methods, we have av- 
eraged the result over all initial impact angles between 
the particle and filament. We find that, for the highly 
ordered magnetic fields required in NRFs, incident elec- 
trons are almost immediately expelled from the filament 
(within seconds) and thus do not significantly contribute 
to the synchrotron emission from these objects. While 
electrons may be able to enter through the edges of a 
filament aligned with the magnetic field, the geometry 
of NRFs suggests that this would occur for only a very 
small (< 1%) fraction of electron impacts, greatly dimin- 
ishing the electron population within the filament. The 
jaggedness of this dataset is statistical in nature, and is 
not relevant for the constraints established here. 

3. SYNCHROTRON EMISSION FROM DARK MATTER 
ANNIHILATION INSIDE OF NON-THERMAL RADIO 
FILAMENTS 

Since magnetic mirroring prevents external electrons 
from propagating into the NRFs, we are motivated to 
study the expected synchrotron contribution of elec- 
trons and positrons which originate from within the 
filaments themselves. The injection of energetic elec- 
trons/positrons through dark matter annihilations is par- 
ticularly interesting for four reasons: (1) the dark matter 
annihilation rate will remain entirely unaffected by the 
peculiar magnetic properties of the NRFs, allowing for 
a rigorous comparison in flux between the filaments and 
the surrounding interstellar medium, (2) the electron in- 
jection spectrum from dark matter will be identical at all 
locations, providing an explanation for the similar syn- 
chrotron spectra observed from multiple NRFs, (3) the 
dark matter density in t he Galactic center is dominated 
by a smooth component ( Diemand et al. 2008 1 , and natu- 
rally explains the intensity observed from all NRFs, and 
(4) a dark matter annihilation scenario at very similar 
energies has already been proposed to explain Galactic 



4 



cente r emission at 7-ray energies ( Hooper & Goodenough 
20111. 



The local dark matter annihilation rate per volume at 
a location x is given by: 



l \ M DM 



(2) 



where p(x) is the dark matter energy density, Mdm is 
the dark matter particle's mass, and (crv) is dark matter's 
annihilation cross section multiplied by the relative veloc- 
ity of the WIMPs. We adopt a dark matter dis tribution 
based on the results of numerica l simulations (Navarro 
et al.|1996| Diemand et al. 2008), and given by p(r) = 
po(r/8.5kpc)~ l:z:> , where r is the distance fr om the 
Galactic center and po = 0.385 Ge V cm -3 (|Catena 



fc Ullio||2010| |Lisanti fc Spergel||2011| . We also adopt 
a dark matter annihilation cross section of (av) — 
3 x 10~ 26 cm 3 s _1 , which is the value predicted for a 
simple thermal relic. 

Assuming approximate cylindrical symmetry for the 
filament geometry, the annihilation rate within a filament 
of length I and diameter w is given by: 



DM 



1.4 x 10 33 s 



33 „-l 




3 x 10- 26 cm 3 /s 
\ / \ 2 



40 pc / I 1 pc 



(3) 



where r is the distance of the filament from the Galac- 
tic center. We note that in many NRFs, the distance 
from the Galactic center changes considerably across the 
length of the filaments - we will discuss the effect of this 
when modeling specific filaments. 

The types and spectra of particles produced in dark 
matter annihilations depend on the details of the parti- 
cle physics model. In order to generate a bright flux of 
synchrotron emission with a spectrum peaking at ~10 
GHz, we will focus on dark matter which annihilates 
dominantly to charged leptons. In particular, we will 
consider a democratic model which annihilates equally 
to e , p ± , and r final states. The nearly instantaneous 
decays of the taus and muons produce lower energy elec- 
trons/positrons, as well as a prompt flux of 7-rays (as 
opposed to 7-rays from the inverse-Compton scattering 
of ene rgetic electrons). It is this pro mpt flux of 7-rays 
which iHooper & Goodenough (2011) find to be consis- 
tent with the excess observed in the Galactic center by 
Fermi-LAT. 

The spectrum of electrons and positrons produced 
through dark matter annihilations withi n a given NRF 
is calculated using the Py thia package ( Sjdstrand et al. 
2001[|Gondolo et al |2004| . In the left frame of Fig.[2j we 
show the injected electron spectrum per dark matter an- 
nihilation for our canonical case of a dark matter particle 
with a mass of 8 GeV, annihilating equally into e + e~, 
p + p~ and t + t~. We note that the electron/positron 
spectrum is very hard, following a spectrum between 
E -0 - 5 and E° between 100 MeV and 8 GeV, although we 
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Fig. 2. — (Left) The spectrum of electrons injected through the 
annihilation of an 8 GeV dark matter particle to e+e~, fi + and 
t + t - (democratically). Note that the contribution from e + e — 
takes the form of a delta function at 8 GeV and is concealed by the 
line denoting the "Total Flux" . (Right) The spectrum of electrons 
from dark matter annihilations before (solid) and after (dashed) 
synchrotron energy losses for an energy loss time of r = 1.0 (as 
defined in Eqn. pi. The direct flux to electrons represents a delta 
function at the 8 GeV mass of the dark matter particle which 
carries approximately 2/3 of the total electron energy. The positron 
flux has an identical energy dependence and lends a factor of two 
to the total lepton flux from dark matter annihilation. 

caution that this spectrum is not a continuous power law. 
The majority (~2/3) of the electron energy is deposited 
in a delta function at 8 GeV following the dark matter an- 
nihilations directly into electrons. These 8 GeV electrons 
will dominate the synchrotron spectrum from NRFs, in 
part due to their shorter synchrotron energy loss time. 
We note that the positron spectrum is identical and lends 
another factor of two to the overall synchrotron flux. 

In order to determine the synchrotron spectrum ex- 
pected from dark matter annihilations, we must also 
model the diffusion of electrons throughout the NRFs. 
Models of the filaments typically require an electron dif- 
fusion timescale similar to the energy loss time of the 



et al. 


1995 


Lang et al. 


1999b) 



This has the effect of 



ening the overall synchrotron spectrum. Through con- 
siderations of the electron gyroradius similar to those 
discussed in § [2j electrons created within the NRFs are 
constrained from effective diffusion perpendicular to the 
ordered magnetic field. In the case of an entirely or- 
dered magnetic field, charged leptons would spiral freely 
along the magnetic field lines until exiting the filament. 
However, in the observed regime containing significant 
ordered and unordered fields, diffusion is expected to be 
significantly more complicated. 

In the case of very low turbulence levels, much work has 
been done within the perturbative framework of quasi- 
linear theory which seeks to calculate the parallel and 
perpendicular diffusion components as a function of the 
power in turbulent modes of the magnetic field with a 
wavenumber resonant with the inverse of the particles 
momentum ( |Jokipii|1966 1971 1. The amplitude of these 
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modes, however, is poorly constrained in galactic simu- 
lations. More recently, numerical simulations have been 
used to analyze the parallel and perpendicular diffusion 
constants in regimes in which th e ordered and uno rdered 
field are co-dominant. Notably, Casse et al. (20021 found 
that in the case of a magnetic field which is approxi- 
mately 80% ordered, the parallel diffusion constant ex- 
ceeds the perpendicular diffusion constant by a factor 
of ^125. Very similar results were lat er obtained for 
the c ase of more energetic cosmic rays (DeMarco et al. 
20071. Since the length travelled by diffusive particles 
can be written as I = \/2Dt, where D is the assumed 
diffusion constant. This implies that perpendicular and 
parallel diffusion will remove particles from the filaments 
on equivalent timescales if the parallel and perpendicu- 
lar extent of the NRFs is approximately 125 to 1. This 
length-scale is similar to the observed length to width ra- 
tio in multiple filamentary structures ( Yusef-Zadeh et al 
2004 1 , implying that diffusion along the length ot the til- 



amentary arcs acts on a similar time frame to diffusion 
across the much smaller NRF width. Since the calcula- 
tion of magnetic field order stands as a lower limit in this 
calculation, it is feasible that perpendicular diffusion is 
in fact entirely irrelevant in the population of NRF struc- 
tures. 

The overall normalization of the diffusion constant de- 
pends sensitively on the length scale o f the turbulen t dis- 
turbances in the magnetic medium ( Jokipii |1966[ ) and 
is highly uncertain. While simulations are able to con- 
strain the mean galactic diffusion constant through ob- 
servations of cosmic ray primary-to-secondary ra tios at 
the solar position (e.g Strong fc Moskalenko|1998 ), these 
simulations do not constrain local diffusion constants, es- 
pecially in magnetically unique regions such as NRFs. 

The synchrotron energy loss time of an electron is given 
by: 



E 
E 



6.6 x 10 



12 



8 GeV \ / 100 fiG 



E 



B 



(4) 



Due to the difficulties of calculating the diffusion con- 
stant within a partially ordered magnetic field, we choose 
a parallel diffusion constant such that electrons remain 
within the NRF for a length of time given by: 



T, 



confinement 



E 
E 



8GcV 



0.33 



E 



6.6 x 10 12 s - 




1.33 



100 fj,G 



(5) 



where r is the ratio of the diffusion timescale for 8 GeV 
electrons compared to their synchrotron energy loss time. 
For example, in the case r = 1.0, 8 GeV electrons diffuse 
out of the NRF on a timescale equal to their synchrotron 
energy loss time. Thus the r parameter can also be seen 
as an indicator of the average synchrotron exhaustion, or 
the average time that an electron generated by dark mat- 
ter annihilation has propagated through the NRF before 
producing the synchrotron emission presently observed. 
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Fig. 3. — The spectrum of synchrotron radiation (in Janskys, 
defined as 10 — 26 watts per square meter per Hz) from electrons 
produced by the annihilation of an 8 GeV dark matter particle 
democratically into leptonic final states in magnetic fields of 50 fiG 
(red, top left), 100 fiG (black, top right), and 200 fj,G (blue, bot- 
tom), with the electron distribution softened during propagation 
times of t = 1.0 (solid), 0.1 (dotted) and 2.0 (dashed). 

The additional factor of E~ - 33 accounts for the e nergy 
dependence of the diffusion constant calculated by |Kol-| 
mogorov (1941). We note that the synchrotron softening 
of an electron spectrum for a given value of r is indepen- 
dent of the magnetic field strength. In Fig. [2 (right), we 
show the spectrum of electrons from dark matter annihi- 
lations after accounting for synchrotron energy losses for 
t = 1.0. 

We are now prepared to calculate the synchrotron spec- 
trum resulting from dark matter annihilations taking 
place within a NRF. In Fig. |3j we plot the synchrotron 
spectrum from dark matter annihilations for magnetic 
field strengths of 50 fiG, 100 fiG, and 200^G and for 
values r =0.1, 1.0 and 2.0. In each case, we predict a 
peak in synchrotron energy at ~1-10 GHz followed by a 
suppression of the synchrotron emissivity at higher fre- 
quencies. In the following section, we will compare this 
prediction to the synchrotron spectrum observed from 
specific NRFs. 

4. COMPARISON TO SPECIFIC FILAMENTS 

In astrophysical interpretations of NRF observations, 
variations in the electron injection spectrum can be in- 
voked to effectively explain the different spectral features 
in each NRF, since the peak of the synchrotron emission 
spectrum depends on the square of the electron energy. 
However, in the case of dark matter annihilations, the 
injected electron spectrum must be uniform in each fil- 
ament. Variations in the observed synchrotron spectra 
may still originate from differences in either the magnetic 
field strength or diffusion timescales of each NRF. These 
effects are relatively weak, however, and would be unable 
to explain extreme variations in the spectral turnover of 
different NRFs. Thus a population survey of the syn- 
chrotron spectra in NRFs remains a powerful diagnostic 
for testing the dark matter interpretation. 

In Table m we have compiled the observed synchrotron 
spectra of the most thoroughly studied NRFs. We find 
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TABLE 1 

Spectral Characteristics of Observed Non-Thermal Radio Filaments 



Name 



Alternative Name 



O: 



4.8GHz 
1.4GHz 



-"4.8GHz 



References 



GO. 08+0.15 
G358.85+0.47 

G359. 1-0.2 
GO.2-0.0 

GO. 16-0.14 
G359.32-0.16 
G359.79+0.17 
G359.85+0.39 
G359.96+0.09 
G359.45-0.040 
G359.54+0.18 
G359.36+0.10 

GO. 15+0.23 

GO. 09-0.09 



Northern Thread 
The Pelican 
The Snake 
Radio Arc 
Arc Filament 

RF-N8 
RF-N10 
Southern Thread 
Sgr C Filament 
Ripple 
RF-C12 
RF-N1 



-0.5 
-0.6 
-1.1 



-0.1 
-0.6 ± 0.1 
0.15 to -1.1* 
-0.5 
-0.5 



-0.5 
-0.8 ± 0.2 
~0.0 
+0.3 
0.0 
-1.0 
-0.9 to -1.3 
-0.6 to -1.5** 



-0.5 to -0.8 
-0.5 to -1.8 
+0.2 to -0.5 



-2.0 
-1.5 ± 0.3 

+0.3 
-0.8 



-0.46 ± 0.32 



Lang ct al. (1999b); LaRosa ct al. (2000) 
Kassim ct al. (iyy9); Lang ct ai. (lyyya) 
Nicholls & Gray (1993); Gray ct al. (iyy5 
|Sofue et j J-iJiJ^p ; |Heich| | [^Un3| l 
|Sofue et al'l {19521) 
LaKosa ct al. (2001) 

|l^aw et al. ( ^Ous flr 

|LaRosa e"t_| i. (2U01); Law ct al. (2008a) 
iLaRosaet al.ljiwOOt 



0.15 



Liszt & Spikcr (iyyt>); Law ct al. (2008a) 
[TJaw et afH2UU8a^ 
|Law et al.| | |2U0gf} 
~w et a l . l ||zuu «a| 
{.eich| | [20D3) | 



Two very different values exist in th e literature for the high frequency spectrum of G359.1-0.2. Gray et al. 1 1995 l cites a value of -0.2 
+ 0.2, while a more recent analysis bylLaw et al. ] {2008b} yields <^^" 8 'q Hz — -1.86 i 0.64 

Spectrum is highly position dependent, but shows a clear trend towards steeper spectral slopes at high frequencies for any given 
position 



the population to be relatively homogeneous, with a hard 
spectrum below ^5 GHz that quickly turns over at higher 
frequencies. The variation in the spectral turnover from 
the hardest NRF (GO.2-0.0, Radio Arc) to the softest 
(GO. 08+0. 15, Northern Thread) is approximately an or- 
der of magnitude. 

In order to test whether magnetic field and diffusion 
timescale variations can explain these spectral and in- 
tensity variations within the highly constrained frame- 
work of a uniform electron injection spectrum, we con- 
sider four NRFs with particularly well measured spec 



tra a nd intensities: GO.2-0.0 (the Radio Arc, Reich 
2003), GO. 16-0. 14 (the Arc Filament , gofue et al. 



1992) 



THTUB+0.15 (the Northern Thread, |Eang°et al.||l999bl 



and G359. 1-0.2 (the Snake, Gray et aT 



1995). TJata were 



extra cted using the Dexter package (Demleitner et al. 
20011, and is shown with the statistical errors and as- 
trophysical background subtraction determined by each 
study. In each case, we calculate the flux at the point of 
maximum 1.4 GHz emission, with the exception of the 
Snake. For that filament, the peak emission corresponds 
to a "kink" in the NRF morphology which shows a spec- 
tral index representative of astrophysical injection. Thus 
for the Snake, the flux is determined at a point where the 
spectral index af ® qh * IS entirely flat, w hich lies at ap- 
proximately 19' in the Gray et al. ( 1995 ) nomenclature. 

In Fig. |4j we provide dark matter fits to the intensity 
and spectrum of these four NRFs. For the Radio Arc 
(top left), we adopt a representative distance from the 
Galactic center of 20 pc, a diameter of 20 pc, a magnetic 
field of 290 /j,G and a diffusion timescale of 7.9 x 10 11 s 
(t = 1.0). For the Arc Filament (bottom left), we adopt 
the same distance of 20 pc, using a filamentary diameter 
of 0.62 pc, a magnetic field of 100 /j,G and a diffusion 
timescale of 1.7 x 10 13 s (r = 2.5). For the Northern 
Thread, we employ a Galactic center distance of 30 pc, a 
diameter of 1 pc, a magnetic field of 50 [iG, and a diffu- 
sion timescale of 2.5 x 10 14 s. Finally, for the Snake we 
adopt a Galctic Center distance of 120 pc, a diameter of 
7 pc, a magnetic field of 100 fiG and a diffusion timescale 
of 1.9 x 10 13 s (t = 2.0). 

One potential mismatch in our best fit parameters con- 



cerns the larger widths necessary for the dark matter 
component to match the intensity of the observed radio 
filaments. While we have produced best fit widths of 7' 
for the Radio Arc, 0.4' for the Northern Thread, 0.25' for 
the Arc Filament, and 2.8' for the Snake, observati onal 
data sup ports smaller widths of 4' for the Rad io Arc (|Re- 
.ich|2003D , 0.07'-0.2' for the Northe rn Thread (jLangeTaT: 
1999b| ), 0.3' for the Arc Filament (ISofue et al ||1992|) and 
0.2'- 1.0' for th e Snake ( Gray et al. 1995 " Anantharama- 
However this mismatch may be expected 



1991) 



iah et al 

for three reasons. First, we have in general attempted to 
match the peak luminosities observed within the NRFs, 
while our simulations calculate the average luminosity 
expected throughout the enhanced magnetic field region 
of the NRF. Changes in the magnetic field intensity and 
structure may confine electrons more effectively in spe- 
cific regions of a NRF, which would then show luminosi- 
ties significantly brighter than those predicted by dark 
matter. Secondly, observations calculating the width of 
the NRFs would be expected to miss the outer regions 
of enhanced magnetic field, where the decreased width 
of the filament would produce only negligible emission. 
Third, the distance to the Galactic center may change 
considerably over the length of the NRF, producing en- 
hanced annihilation rates at regions near the Galactic 
center. For instance, if the Northern Thread is more 
correctly modeled to extend linearly from 20 pc to 45 pc 
from the Galactic center, the total annihilation rate is en- 
hanced by over 50%. However, even these features may 
be insufficient to explain the extremely large diameters 
necessary to explain synchrotron emission in G359.1-0.2 
(the Snake). We note that our analysis does not pre- 
clude astrophysical mechanisms formed within the fila- 
ment from also contributing to the synchrotron emission 
spectrum, which may be the case in this particularly 
bright filament. 

A second issue in this analysis concerns the implica- 
tion that the diffusion constant is significantly lower in 
filaments with particularly soft spectra. These disparities 
are not unique to a dark matter scenario, and have tradi- 
tionally ex isted in astrophysical interpretations of these 
NRFs (cf. |Tsuboi et al^|l995"l |Lang et al.||1999b| ). In 
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both dark matter and astrophysical scenarios, this may 
be understood in a model where the filaments exist as an 
entirely ordered magnetic enhancement superimposed on 
a random diffuse magnetic field of approximately 10 /iG 
which permeates the Galactic center region. The differ- 
ing ratios of the ordered to random magnetic fields (e.g 
80% in the Northern Thread vs. nearly 100% in the Ra- 
dio Arc) would then drive significantly enhanced parallel 
diffusion in the Radio Arc. Alternatively, assumptions 
that the Alfven velocity places an upper limit on the 
speed of electron diffusion implies a diffusion timescale 
which scales as B~ l and would approximately match the 
ratio of diffusion t imescales observed in these two sys- 
tems (Alfven 1942| . We note, however, that this effect 
is not well understood and remains a significant assump- 
tion in our model. Lastly, it is possible that the magnetic 
field structures at the edges of the NRFs are confi gured 
to allow sign ificant reflection of trapped electrons ( |Hey-| 
vaerts et al.||1988[ ). 

Another necessary feature in any dark matter model 
of NRFs concerns the radial dependence of the electron 
injection spectrum. As shown in Eqn. [3j the dark mat- 
ter annihilation rate within a given filament falls off as 
~ r~ 2 - 5 , where r is the distance of the filament to the dy- 
namical center of the galaxy. A quantitative observation 
of the electron injection spectrum in individual filaments 
is difficult, due to the varying lengths, widths, magnetic 
fields, and diffusion constants in the observed filaments. 
However, the distance from the Galactic center to var- 
ious NRFs is thought to span nearly an order of mag- 
nitude, which implies an injection spectrum that varies 
by more than a factor of 300 throughout the NRF pop- 
ulation. This makes the statistical observation of such a 
feature possible, even with extremely crude estimations 
for the astrophysical parameters of individual NRFs. In 
order to examine this necessary trend, we have studied 
the observations of 7 NRFs with i ntegrated fluxes and 



lengths observed at 330 MHz in the LaRosa et al. ( 2000 1 



cat alog, as well as the 13 N RFs observed at 1.4 GHz in 



the |Yusef-Zadeh et al. ( 2004[ ) catalog. In both cases, inte- 
grated fluxes as well as lengths, are provided. We assume 
a constant radial width for all NRF, noting that quoted 
widths for most NRFs fall approximately within a factor 
of two. For this reason, we have removed the Radio Arc 
from our datasets as this assumption is particularly poor 
for that filament. 

The total luminosity of a NRF is expected to depend 
sensitively on its length. In addition to the linear de- 
pendence of the dark matter annihilation rate on the 
length of a filament, longer filaments are expected to re- 
tain electrons for longer periods of time and as a result 
will deposit a greater fraction of their initial energy into 
synchrotron radiation within the filament. In this work, 
we consider three scenarios to account for the influence 
of a NRF's length. First, we we consider the case in 
which electrons are effectively confined and lose their en- 
ergy to synchrotron radiation on timescales much smaller 
than the diffusion timescale (t > 1). In this case the to- 
tal flux in an NRF should depend only linearly on the 
length of the filament. Second, in the case that electrons 
free stream through the filaments on timescales much 
smaller than the synchrotron energy loss time (r <C 1), 
the amount of energy deposited by a single electron into 
the filament is expected to scale with the length of the fil- 
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Fig. 4. — The synchrotron energy spectrum predicted from dark 
matter annihilations (A/qm = 8 GeV, annihilating to e+e — , fi + fi~ 
and t + t - with (<rv) = 3 X 10 — 26 cm 3 s — 1 ) compared to the ob- 
served intensity and spectrum of GO. 2-0.0 (the Radio Arc, top left), 
GO. 08+0. 15 (Northern Thread, top right), G0.16-0.14 (Arc Fila- 
ment, bottom left) and G359.1-0.2 (the Snake, bottom right). The 
magnetic fields, filamentary width, and synchrotron energy loss 
times are shown for the synchrotron match to each filament. 

anient, providing a total flux which scales with the length 
of the filament squared. Finally, in the case that electrons 
diffusively propagate through the filament on a timescale 
smaller than the synchrotron energy loss time (r -C 1 
with D /c -C filament length) the total energy deposited 
by an electron inside the filament will vary as the square 
of the filaments length, providing a total flux which varies 
as the cube of the length of the filament. The cases in 
which the total flux scales with I and I 3 effectively bracket 
the possible degrees of correlation between the length of 
a NRF and it's total flux, while the I 2 case can be con- 
sidered something of a median expectation. 

We first examine the observed dataset at 330 MHz. 
In the left frames of Fig. [5] we plot the flux per unit 
length (top), per unit length squared (middle) and per 
unit length cubed (bottom) as a function of the projected 
distance of each NRF to the Galactic center. In each 
case, we note no significant trend between the distance of 
a given filament from the Galactic center. In other words, 
the distance of a given filament from Galactic center does 
not appear to have significant bearing on its emission 
at 330 MHz, suggesting that astrophysical mechanisms 
(i.e. not dark matter annihilations) are responsible for 
the emission at this frequency. 

The same conclusion is not found at 1.4 GHz, however. 
At this frequency (right), we see a very significant corre- 
lation between the projected distance of a filament to the 
Galactic center and its observed intensity. In particular, 
filaments closer to the Galactic center tend to be consid- 
erably brighter at 1.4 GHz than those farther away. We 
note that for the dark matter halo profile used in this 
paper, we predict a flux which scales with r~ 2 ' 5 , while a 
more generic range of profiles predicts behavior between 
roughly r~ 2 and r~ 3 . 

There are several interesting features of the results 
shown in Fig. [5] First, although the correlation observed 
among the filaments in the 1.4 GHz dataset could have 
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Fig. 5. — Flux per unit length (top), per length squar ed (center), and per length cubed (bottom) for NRFs at 330 MHz (left) |La,Rosa| 
|et al.|2000| ) and 1.4 GHz (right) ( |Yusef-Zadeh et al.|2004| , as a function of the projected distance of each filament from the Galactic center. 
IDrror bars based on the integrated flux, as well as a best-fit linear regression are shown for each frame. Boxed points indicate NRFs which 
are listed in both the 1.4 GHz and 330 MHz datasets. At 330 MHz, there is no clear correlation between the flux and the distance of a 
filament from the Galactic center. At 1.4 GHz, however, those filaments closer to the Galactic center are clearly brighter than those farther 
away. 



plausibly been the result of a selection effect resulting 
from the presence of greater astrophysical backgrounds 
closer to the Galactic center, the 330 Mhz observations 
dispute this reasoning, as does the lack of any bright 
NRFs at 1.4 GHz far from the Galactic center. The dif- 
ferent trends observed in these datasets instead support 
the surprising conclusion that the electron injection spec- 
tra which produces the emission at 1.4 GHz is physically 
distinct from the emission mechanism which dominates 
at 330 MHz. Both trends are naturally explained in the 
case of a light dark matter particle in a magnetic field 
of ~ 100 fiG, as the synchrotron spectrum will peak at 
approximately 1 GHz, with astrophysical processes dom- 
inating the emission at lower frequencies. 

We caution that this relationship is still somewhat 
tenuous, due to the varying NRF widths, magnetic 
fields, three dimensional distances, and diffusion con- 



stants which are not modeled in this analysis, and ad- 
ditional studies will be necessary to better identify and 
refine this apparent trend. If dark matter annihilations 
arc in fact responsible for a sizable fraction of the emis- 
sion observed at 1.4 GHz, improved measurements of 
these parameters should enhance the correlation shown 
in Fig. [31 Furthermore, if astrophysical mechanisms dom- 
inate the low energy synchrotron emission, then we would 
expect the spectral slope in this region to be relatively 
soft. We note that the observed trend is steepest among 
the 1.4 GHz NRFs which were not included in the 330 
MHz analysis, making additional measurements in this 
region necessary to better understand the apparent mis- 
match between low and high frequency observations. Fi- 
nally, we note that no correlation such as that described 
here would be expected in reconnection or shock acceler- 
ation models of NRF. However, models based on a mo- 
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noenergetic electron flux from Sgr A* could plausibly 
lead to a similar relation, although as argued in §|2 such 
electrons are not expected to penetrate into the NRFs. 

5. DISCUSSION AND CONCLUSIONS 

The observed synchrotron emission from non-thermal 
radio filaments (NRFs) in the Inner Milky Way have long 
been difficult to explain with known astrophysical mecha- 
nisms. In this article, we have proposed that dark matter 
annihilations taking place within these filaments could 
produce the nearly monoenergetic electron spectrum nec- 
essary to generate the hard synchrotron emission that is 
observed. In this regard, dark matter annihilations have 
several advantages over proposed astrophysical mecha- 
nisms. First, electrons produced through dark matter an- 
nihilations yield a synchrotron spectrum in good agree- 
ment with the hard spectral index and turnover observed 
from NRFs, without mandating finely tuned magneto- 
hydrodynamic interactions to move electrons indepen- 
dently into equivalent monoenergetic electron distribu- 
tions. Second, WIMPs annihilating into leptonic final 
states, such as those employed at low mass es in order 
to explain the Fer mi-Galactic center excess (Hooper & 
Goodenough|20 11 ) , are predicted to inject electrons into 
the Galactic center region with an energy density very 
similar to the synchrotron signal observed from NRFs. 
Lastly, as shown in Fig. [5j dark matter annihilations 
naturally explain the observed correlation between the 
radial distance of NRFs to the Galactic center and the 
inferred electron injection spectrum in the filaments. 

A dark matter origin of the observed radio emission 
from NRFs yields several concrete and testable predic- 



tions. In particular, the dark matter annihilation rate 
and corresponding flux of injected electrons must show 
approximate sp herical symmetry wit h respect to the 
Galactic center (Diemand et al. 2008), and the injected 
spectrum of electrons must be identical for all filaments 
throughout the Galactic center. This is not the case for 
other proposed astrophysical mechanisms, for which the 
electron injection spectrum can vary from filament to fil- 
ament. Furthermore, we have demonstrated that the in- 
tensity of 1.4 GHz emission is significantly enhanced for 
NRFs near the Galactic center compared to those farther 
away, as would be expected from dark matter distributed 
in a cusped halo profile. Observations which are able to 
independently determine or constrain the magnetic field 
strengths and other characteristics among a population 
of NRFs could be used to further examine any scenario 
requiring a single electron injection spectrum. If dark 
matter annihilations are found not to power NRFs, these 
objects may be used to place stringent constraints on the 
dark matter annihilation rate in the region surrounding 
the Galactic center. New observations of NRFs will thus 
be integral to resolving both the astrophysics and parti- 
cle physics of these unusual regions. 
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